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Abstract 

We discuss the phenomenology of charged Higgs bosons in the MSSM with minimal flavor 
violation. In addition to the constrained MSSM (CMSSM) with universal soft supersym- 
metry breaking mass parameters at the GUT scale, we explore non-universal Higgs mass 
models (NUHM) where this universality condition is relaxed. To identify the allowed pa- 
rameter space regions, we apply constraints from direct searches, low energy observables, 
and cosmology. We find that values of the charged Higgs mass as low as mjj+ ~ 135 GeV 
can be accommodated in the NUHM models, but that several flavor physics observables 
disfavor large contributions, associated with high tan/3, quite independently of MSSM 
scenario. We confront the constrained scenarios with the discovery potentials reported by 
ATLAS and CMS, and flnd that the current exclusion by indirect constraints is similar 
to the expected LHC discovery reach with 30 fb~^ of data. Finally, we evaluate the sen- 
sitivity of the presented discovery potential to the choice of MSSM benchmark scenario. 
This sensitivity is found to be higher in the case of a light {171^+ < rnt) charged Higgs. 

1 Introduction 

Charged Higgs bosons are attractive ingredients in theories which extend the Standard Model 
(SM) Higgs sector with additional fields in non-singlet representations of SU{2)l. In particu- 
lar, the observation of a fundamental charged scalar can provide unambiguous evidence of an 
extended Higgs sector in unfavorable cases when observing a neutral scalar alone would not be 
sufficient. This could be the case e.g. in the Minimal Supersymmetric Standard Model (MSSM), 
with the lightest CP-even Higgs boson having similar properties to the SM Higgs boson, and 
the heavier neutrals escaping detection. 
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Apart from the interest to high energy colhder experiments, the charged Higgs plays an 
important role for several experiments at lower energies. As another mediator of charged cur- 
rent interactions, it can contribute (constructively or destructively) to decay processes without 
requiring flavor structure beyond the CKM framework. This coupling to flavor physics makes 
the charged Higgs particularly useful for constraining indirectly the structure of a Higgs sector 
beyond the SM [1-3]. 

Awaiting the experimental verdict, which hopefully will be delivered by the LHC, some 
questions can be addressed. First of all, which models with a charged Higgs have already been 
probed by existing experiments? Second, what do these experiments imply about the prospects 
of discovering charged Higgs bosons at the LHC? We discuss these questions in the context of 
supersymmetry (SUSY) with a minimal extension of the SM Higgs sector. Earlier work along 
these lines was presented in [4-6]. 

The outline of this paper is as follows. In section [2l the Two Higgs Doublet Model (2HDM) 
of the MSSM is discussed in more detail, as we review briefly the theoretical aspects of Higgs 
bosons in the CP-conserving MSSM with minimal flavor violation (MFV). For the numerical 
studies, we have limited ourselves to two particular classes of MSSM models: the constrained 
MSSM (CMSSM), and models with non-universal Higgs masses (NUHM). These are introduced 
in section [31 

To the CMSSM and NUHM models we apply constraints from direct searches, low energy 
observables, and cosmology in order to identify the parameter regions still open for charged 
Higgs bosons. The constraints are presented in section HJ including numerical results and 
uncertainties which are discussed in some detail. This section extends the recent work [6] by 
including additional constraints, discussing also the heavy charged Higgs, and by allowing larger 
non-universality in the NUHM Higgs mass parameters. 

In section we analyze the impact of applying the constraints by comparing to the collider 
reach for discovering charged Higgs at the LHC. We consider both the case with < 
rrit — rrib, which is of particular interest for early LHC running due to the large number of top 
quarks that will be produced, and mH+ > rrit which becomes interesting at a later stage. Our 
findings are exemplified on simulation results from ATLAS and CMS and we compare to the 
current experimental limit from the Tevatron. We also discuss the dependence of the presented 
experimental reach on the choice of MSSM scenario. Section [HI contains a summary and the 
conclusions. 



with hypercharges F = =pL To avoid flavor changing neutral currents (FCNC) at tree level, a 
discrete symmetry is imposed on the Higgs sector which makes Hi couple exclusively to down 
type fermions, while H2 couples only to up type fermions. This choice corresponds to the type 



2 Charged Higgs Bosons in the MSSM 



To discuss the 2HDM of the MSSM, we introduce the two ^^7(2)^ doublets 
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II 2HDM, which is realized in the MSSM. The Higgs potential then takes the form [7] 



V{Hi, H2) = {ml^ + \fi\^) \Hi\^+ [ml^ + \^Ji\^) lifaP - B^l {e,jH{Hi + h.c.) + 

1 1 (2) 

\9'\H\*Hl\' + \{g^ + g'^){\H,\^-\H,\^)\ 

where e^j is the completely antisymmetric tensor, with eu = 1. The mass parameters m%_^, m^^^ 
may take on negative values in order to break the electroweak (EW) symmetry. A non-zero 
value for Bfi breaks softly the discrete symmetry implemented for the type II model. The 
quartic terms in the MSSM Higgs potential are determined by the EW gauge couplings g and 
g'. In the MSSM, all parameters in the Higgs potential are real, and CV is conserved in the 
Higgs sector at tree level. 

Assigning vacuum expectation values to the neutral Higgs components we define, as usual, 
tan/? = {H2) I (i^i) = V2/V1. Fixing v"^ = vl + to give the correct value for m^, the 
MSSM Higgs sector is completely determined at tree level by two parameters: tan/3 and one 
common mass scale for the Higgs bosons. The physical charged Higgs is defined through 
= —Hf sin 13 + H2 cos (3. Since there is only one free mass parameter, the masses of the 
neutral and charged Higgs bosons are related by tree level mass relations. For our purposes, 
the most important such relation is 

2 2,2 /o\ 

between mH+ and the mass mA of the CP-odd Higgs boson. This relation shows directly the 
quasi- degeneracy of mH+ and tjia for large Higgs masses. 

The charged Higgs phenomenology in the MSSM differs from that of the general 2HDM in a 
few respects: the presence of the EW gauge couplings in Eq. ([2]) ensures that tree level unitarity 
is fulfilled. The additional theoretical constraints lead to mass relations such as Eq. There 
are also new aspects related to the additional states introduced by supersymmetry. Although 
the flavor structure in the MFV framework is left intact, sparticles may contribute indirectly 
to the same flavor observables as charged Higgs through loop effects. New decay chains where 
is produced, or new decay modes for if + into SUSY states, can also appear. 



2.1 SUSY Yukawa corrections 

At loop level, the discrete symmetry is broken in the Yukawa sector, inducing couplings to 
the "wrong" Higgs doublet [8-10]. Some of the induced corrections are enhanced by tan/3, 
and for precision phenomenology it is important to take them into account. Using an effective 
Lagrangian approach, the charged Higgs coupling to fermions is modified as follows [11] 



rr md tan 3 , ^ 

V,,n,,t^P^V-'^^^, (4) 



where m^^ and V^^ refers to experimental quantities. Assuming perturbation theory remains 
valid, |eitan/?| < 1. The correction ii is composed from the two pieces 

= Co + e2(5j3, (5) 
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where the second term, proportional to the top Yukawa couphng yt, only contributes to the 
couplings involving the top quark, i.e. tb, ts, and td. Explicitly, the tan (3 enhanced corrections 
are [12] 



,, = -?^^hJ<^] (6) 

Sit rrin \ mi mi 
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IGvr^ /i y ' |^|2 

where mQ,mij,m'jy are generation dependent soft SUSY breaking scalar masses, and At the 
trilinear stop coupling. The function H2{x,y) is defined as 

IT r \ ^1^^ , y^'^y /ox 

H2{x,y) 



'l-x){x-y) {l-y){y~x) 
For the third generation Yukawa correction, the compact notation 

Ab = 63 tan p (9) 

is sometimes used. Considering the decoupling limit, when the SUSY masses are sent to a 
common high scale, the expression for q becomes particularly simple. Neglecting the second 
part of Eq. ([5]), one obtains q = eo = sign(/i) x a;s/(37r). This limit also gives the simple 
estimate |eo| — 0.01 for the magnitude of the corrections. 

In a renormalization group improved treatment, frib should be renormalized at hr = mH+ to 
include QCD effects to all orders [13], while q is preferentially evaluated at the SUSY breaking 
scale MsusY [10] to avoid large logarithms of the type log(/ii{/MgusY)El 



3 Specification of the MSSM models 

We consider two specific scenarios to illustrate the constraints and collider prospects for charged 
Higgs bosons in the MSSM. Both models assume SUSY breaking mediated by gravity, minimal 
flavor violation (MFV), and conservation of 7^-parity. The first model is the constrained MSSM 
(CMSSM), characterized by the set of parameters {mo, mi/2, Aq, tan/?, sgn(/i)}. Here mo is 
the universal mass of the scalars, mi/2 the universal gaugino mass, Aq the universal trilinear 
coupling, and tan (3 the ratio of the vacuum expectation values of the Higgs doublets, as given 
above. The CMSSM model invokes unification boundary conditions at a very high scale mcuT 
where the universal mass parameters are specified. The masses at the EW scale are determined 
through renormalization group evolution. Additionally, the radiative corrections must generate 
the correct shape of the Higgs potential in order to break the EW symmetry. 

The second model we consider involves non-universal Higgs masses (NUHM). This model 
generalizes the CMSSM, allowing for the GUT scale mass parameters of the Higgs doublets to 

*The recently completed two loop calculation indicates a change of 0(10 — 15%) in £3 [14]. Both the sign 
and magnitude of the correction depends on the MSSM scenario. The renormalization scale dependence is 
significantly reduced compared to the one loop result. 
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Parameter 


min 


max 


note 


mo 


50 


2000 




mi/2 


50 


2000 




Ao 


-2000 


2000 






-2000 


2000 


CMSSM: only sign ± 


rriA 


5 


600 


NUHM only 


tan/? 


1 


60 





Table 1: Parameter ranges used for the CMSSM and NUHM scans. Dimensionful values (all 
parameters except tan/3) are given in GeV. 

have values different from mo, i.e. rriHi 7^ 7^ "^o- These two additional parameters with 
dimension of mass can be traded for two other parameters at a lower scale, conveniently the /i 
parameter and the mass m^ of the CP-odd Higgs boson. The tree level mass relation ([3]) then 
implies that the charged Higgs boson mass can be treated essentially as a free parameter, an 
important difference comparing to the CMSSM. 

To investigate the parameter spaces of the CMSSM and NUHM, we perform scans in a 
random grid generating of order 10^ points for each scenario. The input ranges used for the 
parameters are given in Table [H 

The spectrum of SUSY particle masses and couplings is calculated for each model point 
using SOFTSUSY 2.0.18 [15]. Full 2-loop RGE evolution is employed for all MSSM parameters, 
Yukawa and gauge couplings. The scale at which the universal MSSM parameters are specified 
is determined by unification of the EW gauge couplings. 

In experimental simulation studies, as performed both by ATLAS and CMS, an updated 
version of the m/j-max scenario [16] is used as a common benchmark. This scenario is phe- 
nomenologically defined by the weak scale parameters 

MsusY = 1 TeV 
Ms = 200 GeV 
Mg = 800 GeV 
= 2 TeV 
fi = 200 GeV. 

The gaugino masses have values inspired by gauge coupling unification, and Mi = 5/3M2 tan^ 9iy 
X^^ is defined in the on-shell scheme and relates the amount of mixing in the stop sector to 
the trilinear coupling At through X^^ = At — fi cot p. Compared to the original m/j-max sce- 
nario, the preferred sign of fi has later been changed to positive. As we will see below, this 
accommodates better for the experimental results on the anomalous magnetic moment of the 
muon. 
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4 Constraints 



We use a set of direct and indirect constraints in order to determine the parameter space 
regions allowed for charged Higgs bosons in the CMSSM and NUHM scenarios. Present data 
already provide interesting information on the models, while improvements in the theoretical 
calculations of both the Standard Model and supersymmetric contributions to a number of low 
energy observables increase their predictability. Some recent analyses showing constraints on 
the MSSM parameter space can be found in [4,5,17-21]. 

To obtain constraints on {■mH+, tan jS) we consider a) direct mass limits from LEP and 
the Tevatron, b) flavor data constraints, c) the muon anomalous magnetic moment, and d) 
the dark matter relic density. The flavor observables and [g — 2)^ are calculated with Superlso 
v2.3 [22,23]. A brief description for each observable is given below, and a more detailed account 
of the calculations can be found in [23]. MicrOMEGAs 2.1 [24,25] is used for computing the 
dark matter relic density. For brevity, we illustrate the exclusion by different constraints using 
the NUHM model points only. The combined constraints are presented both for the CMSSM 
and the NUHM models. 

4.1 Direct mass limits 

The non-observation of the charged Higgs boson, or any other SUSY particle, in direct search 
experiments at LEP2 and the Tevatron sets stringent limits on the masses of these particles. 
In Table [2] we have compiled a list of mass constraints from the PDG [26]. Some of these limits 
are subject to auxiliary conditions, e.g. tan/5 < 40, which we take into account consistently. 



Particle 


H+ 


h 


X? xt 


eR 






hi h g 


Mass limit (GeV) 


79.3 


111 


46 94 


73 


94 


81.9 94 


89 95.7 308 



Table 2: Lower limits at 95% C.L. on masses of sparticles and MSSM Higgs bosons. Some 
limits are subject to auxiliary conditions which are not listed in the table. We refer to [26] for 
the complete list and further details on how they are obtained. 

The limit on the mass m/j of the lightest CP-even Higgs boson is very important for con- 
straining also the heavier Higgs bosons, including m//+. This is a consequence of the tree 
level relations between the Higgs masses in the MSSM. In this work we apply the SM limit 
rrih ^114 GeV [27]. Assigning a 3 GeV intrinsic uncertainty in the Higgs mass prediction 
from higher order corrections [28], the value we finally use is rrih > HI GeV. As discussed 
in [19,29,30], it is possible in the MSSM to lower this bound down to rrih — 90 GeV, in partic- 
ular for high tan (3, by reducing the ZZh coupling. The precise value for this limit depends on 
the MSSM scenario, requiring a complete reanalysis of the LEP Higgs search results for each 
model point, which is beyond the scope of the present work. 

Figure [1] shows the effect in the (m//+, tan (3) plane of imposing the direct mass limits on the 
NUHM model points. The scan over the six dimensional parameter space is treated as follows: 
in the left plot, the allowed points are displayed in the foreground. The opposite is true for the 
right plot, where the excluded points are shown in the foreground. In this way the dependence 
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Figure 1: Constraints from direct mass limits on the six dimensional NUHM parameter space, 
projected onto the plane (m//+, tan/5). The left plot shows the allowed points (green) in the 
foreground, whereas the right plot shows instead the excluded points (black) in the foreground. 
The difference illustrates the dependence on the other NUHM parameters. 

on the model parameters which are not shown in the figure can be determined by comparison. 
A large difference between the two plots in a particular region indicates strong dependence on 
the other parameters, whereas a small or no difference indicates weak dependence. 

Some distinct features of the direct mass constraints are visible in Figure [H There appears 
a lower limit of mH+ ^ 135 GeV, independent of tan/5. This value follows directly from the 
mass relation in Eq. ([3]) with rriA = HI GeV. The higher order corrections to this relation 
are typically small in the NUHM models. This is not necessarily true in the more general 
MSSM, where special parameter regions with light squarks (|/i| > AMsusy) allow for large 
mass sphttings in the Higgs sector, even to the degree mH+ < mA — mw [31]. Lowering the 
limit of iTih, as aforementioned, would result in a corresponding shift in the mH+ limit. The 
region with tan /? < 3 is excluded, again as a result of the rrih limit. Except in these two regions, 
the direct mass limits do not provide further constraints on (mj:^+, tan /?) in the NUHM models. 

In the CMSSM, the direct mass limits are not more constraining than in the NUHM models. 
On the contrary, there is a region for high 771^+ and tan (3 which is always allowed. This is 
because in the CMSSM at high tan/3, tuq and mi/2 cannot be too low, and h, Xi and g are 
always sufficiently heavy to avoid the direct constraints. 

4.2 Flavor physics constraints 

Constraints on charged Higgs bosons can be obtained from low energy flavor physics experiments 
by measuring the decay rates of B and K mesons and comparing to the SM predictions. This 
is challenging in many respects: experimentally, because many of the interesting transitions are 
rare, and theoretically, since the predictions often suffer from large hadronic uncertainties. In 
the following, theoretical and experimental uncertainties are added in quadrature. The value 
mt = 172.4 GeV [32] is used throughout. 
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4.2.1 6 ^ S7 

The rare FCNC process b s'j, occurring first at one loop level in the SM, allows for new 
physics contributions from a charged Higgs loop to be of comparable magnitude. Since the 
charged Higgs always contributes positively to this branching ratio, it is an effective tool to 
probe the 2HDM. In the MSSM however, there exist additional contributions from loops with 
charginos and squarks, which may be either negative or positive. Hence the charged Higgs 
constraints in the MSSM are not necessarily as strict as those obtained in the pure 2HDM. 

Following [33] , the theoretical prediction for the inclusive branching ratio of 6 — > 57 can be 
written 



BR(B ^ Xs7) = BR(B ^ X^eP) 



exp 



Vcb 



^[P(Eo) + N(Eo)], (10) 



where P{Eq) and N{Eq) denote, respectively, the perturbative and non-perturbative contribu- 
tions evaluated for a cut E > Eo on the photon energy. C is a semi-leptonic normalization 
factor. The details are described in [33]. The SM prediction for this decay is known to NNLO 
accuracy [33-35] . Using the updated input values of [26] , we obtain 

BR(B ^ Xs7)sM = (3.06 ± 0.22) x 10-^ (11) 

while one would retrieve BR(B — > Xs7)sm = 3.15 x 10~^ with the input values of [33]. The 
combined experimental value of the branching ratio is updated by the Heavy Flavor Averaging 
Group (HFAG) [36]. Their latest result is 

BR(B ^ X,7)exp = (3.52 ± 0.23 ± 0.09) x 10"^ (12) 

To obtain the allowed range for BR(B — > Xs7) in the MSSM, we follow the procedure out- 
lined in [17, 19], in which the intrinsic MSSM uncertainty is added to the uncertainties of the 
experimental value and the SM prediction. The resulting allowed range at 95% C.L. is 

2.15 X 10"^ < BR(B Xs7) < 4.89 x 10"^ (13) 

Another constraining observable which can be extracted from b ^ transitions is the 
degree of isospin asymmetry. It has been shown earlier [17, 37] , that it often provides stricter 
limits on the parameters of different MSSM scenarios than the inclusive branching ratio. The 
isospin asymmetry Aq in the exclusive decay B —>■ K*j is defined as 

_ r(^° ^ x*°7) - r(B^ x*±7) 
°± - r(so ^ K*^^) + r(s± ^ x*±7) ■ ^ ' 

It can be determined from [38] 

Ao = Re(6d-6j, (15) 

where 

K= 'K^'S' f^i^i+^T'^^^l (16) 
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rriH. (GeV) m^. (GeV) 

Figure 2: Constraints from the branching ratio and isospin asymmetry in 6 — > 37 transitions on 
the six dimensional NUHM parameter space, projected onto the plane (mj:/+, tan/3). The left 
plot shows the allowed points (green) in the foreground, whereas the right plot shows instead 
the excluded points (red) in the foreground. 

are spectator quark dependent terms. We refer to [37, 38] for the definition of the different 
terms, and their expressions in terms of Wilson coefficients. 

Combining the most recent experimental values of BaBar [39] and the results of Belle [40], 
including the experimental and theoretical uncertainties, the allowed range 

- 1.7 X 10"2 < Ao < 8.9 X IQ-^ (17) 

is obtained at 95% C.L. [23]. 

In the same way as for the direct constraints, the results of applying the 6 — 57 constraints 
to the NUHM model points are displayed in two separate plots in Figure O The figure includes 
both the constraints from the branching ratio and the isospin asymmetry combined. Low values 
for mH+ and tan j3 simultaneously are excluded, regardless of the other NUHM parameters. As 
expected however, it is possible in most of the parameter space to balance out the contributions 
from charged Higgs and from chargino-squarks against each other to a sufficiently high degree 
to be consistent with the SM. Observables based on b s'j transitions thus represent high 
sensitivity to the MSSM parameters. 

We find that it is easier in NUHM models to avoid the constraints from b s'-f than it is 
in the CMSSM. The isospin asymmetry is more restrictive in the CMSSM than the constraint 
from the branching ratio, and the main allowed region in the CMSSM is obtained for high 
tan p. 

4.2.2 Bu TUr 

In contrast to the b ^ sj transitions, where the charged Higgs participates in loop diagrams, 
the process Bu — > can be mediated by already at tree level. Since this decay is helicity 
suppressed in the SM, whereas there is no such suppression for the scalar exchange in the 
limit of high tan/3, these two contributions can be of similar magnitude [41]. 
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rriH. (GeV) m^. (GeV) 

Figure 3: Constraints from BR [B —>■ rz/^) on the six dimensional NUHM parameter space, 
projected onto the plane {mH+, tan (3). The left plot shows the allowed points (green) in the 
foreground, whereas the right plot shows instead the excluded points (blue) in the foreground. 



The leading order SM prediction for this decay is 



BR{B^ ^ Tu,)sM = "^"1 msmi 1 - -f , (18) 



while the new physics contribution from is expressed through the ratio [42] 



nMSSM _ BR(i?u T^'r)MSSM 



"TUt 



BR(5„ -> rz/^)sM 



y 1 + eo tan (3 



H 2 



(19) 



The leading SUSY-QCD corrections are included in this expression through eo- Using fs = 
200 ±20 MeV [43], and the combined value |Kfe| = (3.95 ±0.35) x 10"=^ [26], the SM branching 
ratio evaluates numerically to 

BR(5„ ^ rz/^)sM = (1-10 ± 0.29) x 10"^ (20) 

The SM prediction is compared to the current HFAG value [36] 

BR(B„ ^ rz/^)exp = (1-41 ± 0.43) x 10"^ (21) 

by forming the ratio 

R7Z - ^ """^V" = 1-28 ± 0.38. (22) 

This ratio still suffers large uncertainties from the determination of \Vub\i since different mea- 
surements of this quantity are incompatible. The constraints obtained should therefore be 
treated merely as an indication, rather than as a strict limit on the same level as, for example, 
6 — *■ S7 transitions. Requiring i?^^^^ to be within 2 a of R^Jl provides the following allowed 
range: 

0.53 < ^ < 2.03 (23) 
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itIh- (GeV) rriH. (GeV) m„ (GeV) 



Figure 4: Constraints on NUHM parameters obtained from BR {B^ — > rur) for three values of 
\Vub\ from [26]. Alfowed points shown in green and excluded points in blue. From left to right 
\V^l^^\ = (4.12 ±0.43) X 10-3, \V^r^\ = (3.95 ±0.35) x lO^^, and \VX^\ = (3.5 ± 0.6) x 10"=^. 



where we have estimated the residual MSSM errors to be negligible. 

The resulting constraints from Bu —>■ are illustrated in Figure [3l It can be seen that a 
fairly large share of the available parameter space is affected. The allowed points fall in two 
disjoint regions. At low mH+ and large tan /5, the contributes twice the SM amplitude with 
opposite sign [see Eq. f|T9l) ]. The exclusion power varies only weakly between different MSSM 
models, as shown by the similarity of the left and right plots in Figure [31 This results from the 
fact that the only source of significant MSSM scenario dependence in this tree level observable 
is through the eo corrections. That this is a tree level observable also means that the results 
can be carried over essentially unchanged to the CMSSM, or to any MSSM model with MFV 
and 7?.-parity conservation. 

Treating the difference in \Vub\ determination as a theoretical uncertainty would eliminate 
the constraint from BR (i?„ tUt). As an alternative, we present in Figure H] the resulting 
constraints on (mj|/+, tan/3) obtained using three different values for \ Vub\, corresponding to two 
separate determinations and to the combined value used for Figure [3] above. 



4.2.3 B Drur 



Compared to Bu rur, the semi-leptonic decays B —>■ Diu [44-46] have the advantage of 
depending on \Vcb\, which is known to greater precision than \Vub\- In addition, the BR(i? 
DrUr) is about 50 times larger than BR(i?u rur) in the SM. The experimental determination 
remains however very complex due to the presence of at least two neutrinos in the final state. 
The branching ratio, including the SM and charged Higgs contributions, can be obtained from 
[46] 

dT{B Diu) Gl\V^b?m% 



dw 



1927r3 



-pv[w) 



X 



1 - 



m 



B 



tiw) 



nib 



tan^/3 



{nib — m'c)'mH+ 1 ± tan (3 



psiw) 



(24) 



where the kinematic variable w = vd ■ vb is written in terms of the meson four velocities, and 



tiw) 



m 



B 



± m|) — 2w mDmB- The definitions of the scalar and vector form factors {ps and 
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Figure 5: Constraints from BR [B DtUt-) on the six dimensional NUHM parameter space, 
projected onto the plane (m/^+, tan/3). The left plot shows the allowed points (green) in the 
foreground, whereas the right plot shows instead the excluded points (orange) in the foreground. 

Pv, respectively) can be found in [46]. To reduce some of the theoretical uncertainties, the ratio 

_ BR (B ^ Dru,) 
^""^ = BR (B Deu^) ^^^^ 

is considered, which is expected to be sensitive to charged Higgs contributions through the final 
state with a r lepton. The SM prediction for this ratio is 

eg^, = (29±3)xl0-^ (26) 

where the main uncertainty comes from the form factors [46]. The most recent experimental 
result by the BaBar collaboration is [47] 

= (41.6 ± 11.7 ± 5.2) X 10-2. (27) 

To derive the allowed range for this observable in the MSSM, the theoretical and experimental 
results are combined. Including also the enhancement of the form factor uncertainties by the 
presence of the charged Higgs contribution in Eq. (^^, we use the following interval at 95% 
CL. in our analysis 

15.1 X 10'^ < ^Deu < 68.1 X 10"^ (28) 

The results of applying the B DrUr constraint to the NUHM points are shown in 
Figure [5l A narrow strip at large tan j3 is excluded for mH+ < 200 GeV. This excluded region 
provides excellent complementarity to the Bu — > TUr constraints obtained above, as most of 
the parameter space region still allowed at low mH+ is covered. Similar (small) MSSM model 
dependence through eo applies in this case as for Bu — > tUt, and the result is unchanged in the 
CMSSM. 
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4.2.4 



As a final B meson decay we investigate the rare process Bg yu^/i^, which has so far not been 
observed experimentally. At high tan jS, the MSSM contribution to this process is dominated 
by the exchange of neutral Higgs bosons. We therefore expect indirect constraints on mH+ and 
tan/3 from the MSSM mass relations. The BR(i?s — > fi~^fi~) can be expressed as [48] 



647r3 



fBjBMllVtbV, 



* |2 



, ^ _ Ami 



X 



Ami 
Ml 



Mi\Cs\ 



+ 



Ml 



CpMb, 



2 0,"^^ 



M 



Bs 



(29) 



where the coefficients C5, Cp, and Ca parametrize different contributions. Within the SM, Cs 
and Cp are small, whereas the main contribution entering through Ca is helicity suppressed. 
In the MSSM, both Cs and Cp can receive large contributions from scalar exchange. The Bg 
decay constant /^^ = 245 ± 25 MeV [43] constitutes the main source of uncertainty in this 
expression. The SM prediction is 



BR(5, 



)sM = (3.2±0.5) X 10- 



while the current experimental limit, derived by the CDF collaboration, is [49]: 

BR(5, ^ ^+/i") < 5.8 X 10"^ 



(30) 



(31) 



at 95% CL. The experimental limit is thus still an order of magnitude away from the SM 
prediction, allowing for substantial SUSY contributions. Including theoretical uncertainties, 
we compare the MSSM prediction to the upper limit at 95% C.L. 



BR(5, ^ /i+/i~) < 6.6 X 10" 



(32) 



The resulting constraints from Bg ^^^^ are shown in Figure [6l The indirect dependence 
of this constraint on mH+ is clearly visible; in the right plot the possible exclusion is seen to be 
quite effective, whereas there is a large transition to allowed points going to the left plot. Hence 
the dependence on the MSSM scenario is large, and the constraints on {m tan f3) become 
dependent on the masses of the sparticles, for example the charginos. The only region which is 
almost completely excluded is for very small mH+ and large tan (3. 

In the CMSSM, Bg —>■ fj^^fi^ removes most points with tan/5 > 50, complementing the 
direct constraints and b — * 57. 



4.2.5 K 



Another decay process, which has many similarities to Bu — > tUt, is K ^ //i/^. This decay can 
also be mediated by the charged Higgs at tree level, although in this case the contribution is 
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Figure 6: Constraints from BR(i?5 /^^A* ) on the six dimensional NUHM parameter space, 
projected onto the plane (m//+, tan/3). The left plot shows the allowed points (green) in the 
foreground, whereas the right plot shows instead the excluded points (yellow) in the foreground. 



reduced by the coupling of to lighter quarks. In order to reduce the theoretical uncertainties 
from Jk, the ratio of partial widths 



r(vr /iz/^) 



ud 



K 



X 



1 - 



m 



Ml. 



tan^ /3 



1 + eo tan /? 



;i + 5c 



(33) 



is usually considered. Here 5em = 0.0070 ± 0.0035 is a long distance electromagnetic correction 
factor. As suggested in [50], we study instead the quantity 



R 



£23 



0^ 



£3 



Vudi'^£2) 



(34) 



Here VusiKa) refers to Vus as measured in leptonic decay of K with i particles in the final state 
(two leptons and a number of pions), and similarly for Vud- The O"*" — > O"*" denotes nuclear beta 
decay. In the SM i?^23 = 1? while the contribution from charged Higgs in the MSSM attains 
the simple form 



R, 



£23 



m 



Ml. 



tan^ P 



(35) 

1 + eo tan /? ■ ^ ^ 

Using md/rris = 1/20 [26], the MSSM prediction can be directly compared to the experimental 
value [50] 

Rm = 1.004 ± 0.007. (36) 

In the extraction of this value, the ratio fx/fn has been fixed to the value fx/ fn = 1.189±0.007 
obtained from lattice QCD using staggered quarks [51]. It should be noted that the uncertainty 
thus obtained for i?^23 is most probably overly optimistic. Indeed, many approaches exist to 
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Figure 7: Constraints from K — * /iz/^ transitions on the six dimensional NUHM parameter 
space, projected onto the plane (m//+, tan/5). The left plot shows the allowed points (green) 
in the foreground, whereas the right plot shows instead the excluded points (gray) in the 
foreground. The value fx/ fw = 1-189 ± 0.007 was used in obtaining these constraints. 



determine fx/fn, and some reservation remains about staggered fermions [52]. If, for example, 
the value fx/ fn = 1-205 ± 0.018 (obtained using the domain wall formulation [53]) is used 
instead, the BR [K — > fiu) provides no constraints on the studied NUHM parameters. We 
therefore stress that the constraints obtained from fl36|) should serve only as an indication. 

As can be seen from Figure [71 more precise estimates of this observable would be very 
useful in constraining the region with low mjy+ and large tan/3. Such results also provide 
complementarity to the different B decays discussed above. The MSSM model dependence of 
the obtained limit is weak, as expected for a tree level observable. A similar exclusion region 
is obtained in the CMSSM. 



4.3 Muon g-2 

The anomalous magnetic moment of the muon = [g — 2)/2 receives non-zero contributions 
from radiative corrections. It has been determined to high precision both theoretically and 
experimentally, and can therefore be used to probe new physics effects, including the MSSM. 
The latest measured value for based on e+e~ data is [54] 

a^fP = (11 659 208.0 ± 6.3) x 10"^°. (37) 

Comparing this precise measurement to the SM prediction [55] 

aj^ = (11 659 178.5 ±6.1) X 10"^° (38) 

leads to the discrepancy 

Sttf, = a;fP - af^ = (29.5 ± 8.8) x 10"^°, (39) 



15 





100 



200 300 400 
rriH. (GeV) 



500 600 




100 200 



300 400 
rriH. (GeV) 



500 600 



Figure 8: Constraints from the muon anomalous magnetic moment 5a^ on the six dimensional 
NUHM parameter space, projected onto the plane (mj:/+, tan (3). The left plot shows the allowed 
points (green) in the foreground, whereas the right plot shows instead the excluded points 
(purple) in the foreground. 

corresponding to a 3.4 cr deviation from the SM. The 95% CL. allowed range, including uncer- 
tainties from two loop SUSY corrections which have not been included, is: 

11.5 X 10"^° < 5a^ < 47.5 x 10"^°. (40) 

In the MSSM, the discrepancy in 5a^ can be accounted for by the contributions from loops 
with exchange of neutralinos-smuons and charginos-sneutrinos. The result of imposing the 
constraint fHOj) to the NUHM points is presented in Figure [8], which shows no correlation 
between and (mj|/+, tan/3). For the CMSSM, the constraint acts differently, and such 
a correlation is observed. 

It is well-known that the sign of the MSSM contributions to 5a^ is directly coupled to the 
sign of the /i parameter. We illustrate this fact for the NUHM models in Figure [H In the 
following, we will therefore take the constraint (HOj) as a requirement of positive /i values. For 
a recent and more detailed discussion of this issue, see [56] . 

Since the sign of the leading contribution to A;, is determined by /i [compare Eq. ([2])], the 5a^ 
results can also be used to constrain the favored range for these corrections. This is illustrated 
for the NUHM points in Figure [91 If the constraint yU > is taken as an a priori requirement 
on the MSSM model, it means a reduced model sensitivity for all charged Higgs observables at 
high tan /5. 

4.4 Dark Matter density 

When neutral, the lightest supersymmetric particle (LSP) is a suitable candidate for the cold 
dark matter of the universe. In the MSSM, the dark matter is therefore usually expected 
to consist of the lightest neutralino. The latest 5-year WMAP data [57] provides a precise 
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Figure 9: Dependence on the fi parameter of Sa^ (left). Correlations between contribution to 
Sa^ and A;, (right). 

experimental determination of the dark matter density: 

rij^^h^ = 0.1143 ± 0.0034, (41) 

which can be compared to the theoretical calculation of the LSP relic density. Assigning a 
residual 10% theoretical error [58] in the prediction, we obtain the allowed interval at 95% CL. 

0.094 < ^]DM/i^ < 0.135. (42) 

However, as was shown in [59], it is not safe to use the lower bound in the above interval due 
to cosmological uncertainties in the era prior to Big Bang Nucleosynthesis. Also, dark matter 
can be composed of different components in addition to the LSP, which would again falsify the 
lower limit. We therefore discard the lower bound and use only 

^Duh^ < 0.135 (43) 

to extract constraints from the relic density. 

Requiring a neutral LSP, and that the constraint ( H3i) is satisfied. Figure fTOl shows the results 
on the {mH+, tan f3) plane. The conclusion we draw from this figure is that the cosmological 
constraints on the rehc density do not lead to distinct constrained ranges in (m//+, tan/3). This 
is also the case in the CMSSM. 

4.5 Combined constraints and limits 

In Figure [TT] we show a combination of constraints applied to the NUHM model points. The 
result is projected on the (m//+, tan (3) plane. The constraints are applied in the order indicated 
by the legend, and the first constraint by which a certain point is excluded determines its color. 
Points which are not excluded by any constraint are termed allowed and displayed in the 
foreground to indicate the parameter regions still open for H~^. Only points with /i > and 
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Figure 10: Constraints from the density of dark matter on the six dimensional NUHM parameter 
space, projected onto the plane (mjy+, tan/?). The left plot shows the allowed points (green) 
in the foreground, whereas the right plot shows instead the excluded points (light grey) in the 
foreground. 

a neutral LSP are shown. From the figure, we note that the allowed points fall in a distinct 
region, forming a triangular shape in the lower half plane. The region of allowed points shares 
a diffuse boundary with that excluded by Bu rur transitions. This diffuseness is the result 
of eo variations. 

Taking these constraints into account, we see that charged Higgs masses down to mH+ — 
135 GeV can be accommodated, with the lowest masses allowed for intermediate tan/5 ~ 
7-15. For higher tan/5, the combined constraints follow the exclusion by tUt-. The 

combined results are therefore particularly sensitive to the uncertainties associated with this 
decay channel. The allowed region is given approximately by tan (3 < (m/f+/10 GeV). Above 
the large region excluded by B^ — »■ rz/^ appears a smaller, mostly yellow, region containing 
points excluded by one or more of the constraints from Bs f^~^fJ'~, B DrUr, and K — ^ /xz/^. 
Although the general indication is that is excluded in this region, a more conservative 
conclusion would be that it is still possible to find points excluded by only one of the three 
constraints. Since these can each be questioned on different grounds, there might be some room 
open for alternative interpretations. 

As a final remark on the NUHM model, we note that the constraints from the well- 
established 6 — > S7 transitions are not particularly strong in this class of models. The primary 
exclusion region for 6 — > 57 at low mH+ has already been excluded by the direct constraints, and 
the remainder of the excluded points do not form a distinct exclusion region in the {■mH+, tan (3) 
plane. 

Turning now to the CMSSM, we show in Figure [12] the same combination of constraints 
applied to this model. In the CMSSM, the mass scale for the Higgs bosons is not a free param- 
eter, but it is fixed by the universality assumptions at the GUT scale through RGE running. 
Lower masses for the heavy MSSM Higgs bosons requires the tan j3 enhanced contribution from 
l/f, to cancel the always large RGE effects from the top Yukawa coupling yt. The result of 
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(GeV) 

Figure 11: Combined exclusion in NUHM models by different constraints, as described in the 
text. The constraints are applied in the order they appear in the legend, and the color coding 
corresponds to the first constraint by which a point is excluded. All points have > and a 
neutral LSP. 




itIh. (GeV) 

Figure 12: Combined exclusion in CMSSM models by different constraints, as described in the 
text. The constraints are applied in the order they appear in the legend, and the color coding 
corresponds to the first constraint by which a point is excluded. All points have /i > and a 
neutral LSP. 
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this balancing is seen in the figure, where the distribution of points reveals a clear preference 
for large mH+, and where smaller values for 'mH+ are only obtained in combination with high 

tan/5 > 50. 

The combined constraints in the CMSSM work similarly as for the NUHM in limiting the 
parameter space available for the charged Higgs, with some important differences. The direct 
limits are very effective in ruling out parameter space regions with low-intermediate tan /3 and 
mH+ < 400 GeV. Constraints from Bg — > H'^H" exclude high tan/3 up to mH+ — 500 GeV. 
In the CMSSM in general |Ab| < 0.5 which reduces the MSSM model dependence. The same 
caveats discussed above apply when interpreting the constraints in the CMSSM for low mH+ 
and high tan/3. However, even with some of the more uncertain constraints removed, the 
allowed region in {mH+,ta.n(3) still appears to be minimal. Nevertheless, this region around 
mH+ — 200 GeV is very interesting for early LHC running. 



5 Charged Higgs at hadron colliders 

There are two primary modes through which charged Higgs bosons can be produced at hadron 
colliders. When the charged Higgs is light < rrit — mb), it is kincmatically accessible in 

the decay t bH^ of the top quark. The decay width in the large tan (3 limit reads [10] 

In favorable cases with small mH+ and large tan /3, the 2HDM branching ratio (A;, = 0) may 
reach values up to BR(t — > bH^) ^ 0.3-0.4. The SUSY corrections entering through can 
modify the pure 2HDM value substantially. 

When the charged Higgs becomes too heavy to appear in the decay of on-shell top quarks, 
the favored mode of H'^ production is in association with a single top quark. For this mode of 
production, proper matching is required between the twin processes gg/bg [60]. The full NLO 
calculation is also available [61-63], and we use a parametrization of this cross section for the 
comparison to experimental results below. The 2HDM cross section is then augmented with 
the appropriate tan/3 enhanced corrections proportional to 1/(1 + A^)^. 

The decay of the hght if"*" proceeds mainly through one of the two channels if + — > t'^Vt 
which dominates for tan/3 > 2, or — > cs, becoming important for smaller tan/3 values. 
The tan/3 enhanced corrections to — > r'^v^ are negligible, since there are no SUSY-QCD 
corrections to the leptonic final state. Consequently, the width is given simply by 

T{H+ ^ r+ur) = -l-^mn^ml tan^ /3. (45) 

For a heavier charged Higgs, the decay — > ih opens up, and quickly overtakes — > r'^Vr 
as the dominant mode. In the large tan /3 regime, the partial width 

TV (i^Wth?'^c ( ITT? \^m?tan^/3r .1 
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Figure 13: Total branching ratio for into SUSY particles for the NUHM points, projected on 
the {171^+ ,taii (3) plane (left). The branching ratio multiplied by the charged Higgs production 
cross section (right). 

is proportional to m^, and thus affected by SUSY corrections in the same way as r(t hH^). 
For all numerical evaluations of the branching ratios we use HDECAY [64], which includes 
both QCD and MSSM corrections in a consistent fashion. 

5.1 SUSY decay modes 

When allowed by the kinematics, the charged Higgs may decay to SUSY partners of the SM 
particles. Figure [13] shows the total branching ratio for all SUSY decays of in the NUHM 
models. We see that the largest branching ratios are obtained in the intermediate region 
tan/5 = \Jmt/mb ~ 7, where T[H^ tb) has a minimum, and the detection of charged Higgs 
through the standard decay channels is most difficult. On the other hand, in the right panel of 
Figure [12] we show the branching ratio multiplied with the production cross section. This 
result comes out independent of tan (3. 

We expect the main SUSY decay modes to be xtx'j [65]. The experimental signatures 

for these modes depend on the further decay of the sfermions, where leptonic decays are likely 
to be required to suppress an otherwise overwhelming QCD background. For direct decay into 
the lightest chargino and the LSP, the final state may contain a single charged lepton and 
large missing from Xi ~^ Xi'^'^^e- In the case where the charged Higgs decays to heavier 
charginos or neutralinos, the final state can become more involved. However, a promising 
generic signature is that based on three charged leptons and missing px [66]. 

For the allowed points in the NUHM models we have verified that the chargino-neutralino 
decay modes are completely dominating, and that remaining SUSY decays can be neglected. 
Figure [H] shows the sum of a x BR {H~^ X^X^) for all chargino-neutralino channels to- 
gether, with the constraints of section H] applied. We observe that in particular the direct mass 
constraints, and the constraints from b 57, rule out points with low m o and m+, which 
kinematically would give the highest number of events in the SUSY decay channels. Keeping 
in mind that Figure [TJ] shows only the sum of all chargino-neutralino channels, and that no 
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Figure 14: Charged Higgs production cross section times BR (if + XiXf)^ summed over all 
neutralino and chargino species. Direct and indirect constraint are applied as before. The color 
coding in the left plot agrees with that in Figure 11. In the right plot, only points allowed 
by the constraints are shown, and the color corresponds to ex x BR (in fb). Both plots have 
requirements of > and a neutral LSP. 

branching ratio into one lepton or three lepton final states has been applied, the total cross 
section is of the order of a few 100 fb for the most promising parameter space points allowed 
by the constraints. 

5.2 Tevatron results 

At the Tevatron, CDF [67, 68] and D0 [69] experiments have searched for light charged Higgs 
bosons in the decay of top quarks. The searches have been performed both in the —>■ rVj. 
and — »• cs channels , where the former is of course more interesting in the MSSM. The 
current best limit in the r^Vr channel is obtained by D0 [69] using 1 fb~^ of data. We 

show the model independent limit on BR {t H^b) from this search in Figure [151 assuming 
T^Vr saturates the full width of . 
Recently, there has also emerged D0 results on a search for heavy charged Higgs in the 
— > th channel [70], but with limited sensitivity to the 2HDM (II) at this point. 

5.3 LHC prospects 

The kinematic range of the LHC will allow experiments to search both for light and heavy 
charged Higgs bosons. As discussed above, a heavy charged Higgs would preferentially decay 
through th. However, this channel has proven experimentally challenging. The decay 

mode of primary interest is therefore t^v^ also when mH+ > rrit, even though typically 

the BR{H+ r+z/^) = 10-15% in the hmit of high m^^+. 

To determine the prospects for the LHC experiments to discover the charged Higgs boson 
in the MSSM models under study, we confront our model points with the experimental reach 
for a 5(7 discovery obtained by ATLAS [71] and CMS [72,73] through simulations. For both 
experiments, a full detector simulation is used, and systematic uncertainties are included. The 
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discovery reach is reported for an integrated luminosity of 30 fb^^, corresponding to three years 
of LHC running at "low luminosity" . 

The result for the NUHM model points is presented in Figure [151 showing the light and 
heavy 'mH+ cases separately. For the light mH+, the bulk of the NUHM points are accessible 
already with 30 fb~^. Since the points favored by the constraints lie close to the kinematic limit 
mH+ ^ rrit — m^, dedicated studies are required to carefully evaluate the discovery prospects 
in this mass region. Figure [15] reveals that the models excluded by the different constraints 
are also those which have the highest cross section and branching ratio for charged Higgs 
production. This is simply a result of the universal dependence on tan^ (3 / Tn?^+ shared by most 
observables in the high tan/3 limit. 

Having compared the NUHM models to the experimental discovery reach in a model inde- 
pendent way, we now consider the interpretation in the (mj:/+, tan/3) plane. The experimental 
results [71, 74] are presented in the m/j-max scenario, described in section [31 A comparison 
between the experimental results and the NUHM points is given in Figure [161 The ATLAS 
results in this figure are obtained from a combined discovery contour, whereas the CMS results 
are reported as two contours for light and heavy separately. For ATLAS we also include a 
projected exclusion limit reported at the 2 a level [71]. Figure [T6l illustrates even more explicitly 
than Figure [151 the correspondence between the region with highest discovery reach for collider 
experiments and the most powerful exclusion by indirect constraints. As a side remark, it 
should be mentioned that mH+ < 123 GeV is already excluded in the m/j-max scenario by the 
direct limit on irih [29]. The use of m^-max as a benchmark scenario in this region is therefore 
somewhat questionable. 



5.4 MSSM model dependence 

As a final point we discuss the sensitivity of the presented experimental results to the choice 
of MSSM benchmark scenario. This issue was recently discussed for CMS [74] in the context 
of m/j-max scenarios with different choices for the parameter. For any sub-dominant decay 
channel of a heavy charged Higgs, such as t^t-'t, the effects of the bottom Yukawa 

corrections cancel to a large extent between the production and the decay. To see this, we 
consider the corrected cross section times the branching ratio 

. X BR (H^ -.r^.,)^ A, . ,4r) 

where ctq is the cross section obtained in the pure 2HDM, r^;, the uncorrected width for tb, 
and r the uncorrected total width. Fx refers to any decay mode which is not r+z/,- or tb. 
Equation (HTj) shows that, when F,- and Fx are both small with respect to Ttb, the combined 
At) correction in this channel is second order in BR(if^ —>■ r^Ur), therefore typically less than 
10-15% even for large values of |Afe|. 

To assess the model sensitivity of the ATLAS results presented in [71], we evaluate whether 
the NUHM models could lead to a 5 cr discovery of the charged Higgs boson. For each model 
point, we determine either BR(t — »• bH~^) or cr(pp tH^) as appropriate for the value of 
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Figure 15: Model-independent experimental discovery contours interleaved with NUHM model 
points. Color coding corresponds to exclusion by different constraints (see legend) which are 
explained in the text. Points allowed by the constraints (green) are always displayed in the 
foreground. A neutral LSP and > is required. 
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Figure 16: Experimental discovery contours for the m/i-max scenario, interleaved with NUHM 
model points projected on the (mi^+,tan/3) plane. Color coding corresponds to exclusion 
by different constraints (see legend) which are explained in the text. Points allowed by the 
constraints (green) are drawn in the foreground. A neutral LSP and /z > is required. In the 
low mass case no ATLAS contour at 95% CL is visible, since the reach for exclusion covers the 
whole plane. 
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Figure 17: Sensitivity of ATLAS discovery potential to A;, corrections. Points correspond to 
NUHM models which can be discovered (black), which can be discovered only because of A;, < 
corrections (blue), and points which cannot be discovered because of Af, > corrections (red). 
Both positive and negative values for /i are considered. 

ntH+i followed by the BR(iJ+ t+z/^). This is done both with the MSSM Af, corrections 
applied, and for a fixed A^ = corresponding to a pure 2HDM (II) with the same values for 
(mj^+, tan/3). We then compare the two numbers obtained to what is required for charged 
Higgs discovery with 30 fb~^. In Figure [T7| we show the result of this comparison. Points for 
which the cross section (branching ratio) is large enough for a 5 a discovery with the standard 
Aj, corrections included are shown in blackEl The subset of black points which would not be 
discovered with Af, = are shown in blue. Finally, points which are hidden from discovery 
because of the A^ correction are shown in red. The red points are such that they would be 
accessible with A;, = 0. 

As seen from the left plot in Figure [T71 the BR {t hH^) can be altered quite significantly 
by Ab corrections, resulting in a pronounced MSSM model dependence. There exist (red) points 
for such high values as tan /5 = 50 which do not allow a 5 a charged Higgs discovery. These 
points have large values of A5 ~ 1, thus correspond to large and positive /i. In the NUHM 
models at high tan (3, the resulting distribution of theoretically allowed points is not uniform in 
/i, but has a bias towards positive values. This model effect explains the dominance in number 
of red over blue points in Figure [T7i 

For a heavy charged Higgs the conclusion is different, as seen from the right plot of Figure [T71 
As expected from Eq. P7|) . the dependence on A;, is much milder for this case, resulting in a 
more concentrated distribution of red points. The experimental discovery contour is therefore 
fairly stable with respect to A^ variations. The absence of blue points at high tan/3 again 
results from the model preference for positive in this region. 

*Since we extrapolate the discovery contour in Figure fTSI above ^1^+ = 150 GeV, the distribution in Figure fTTI 
of points allowing a 5 cr discovery does not exactly follow the ATLAS contour of Figure [16] in this mass region. 
This visual difference is of no importance to our conclusions. 
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6 Summary and Conclusions 



Charged Higgs bosons are of special interest since they can provide definite signatures for physics 
beyond the Standard Model. In this paper we have analyzed the charged Higgs phenomenology 
in the constrained MSSM and in models with non-universal Higgs masses. 

To investigate the experimental prospects for charged Higgs discovery, we have first exam- 
ined direct and indirect constraints from a fairly complete set of fiavor physics observables, 
calculated with the publicly available program Superlso [23]. We used results from b 57, 
Bu T^T, B — > DtUt, Bg A*"*"/^", and K — > iiv^ transitions, together with the muon (^f — 2) 
and cosmological constraints. In this manner we have identified the allowed regions for mH+ 
and tan/9 at 95% C.L. The combined constraints from the fiavor observables exclude the region 
with low mu+ and large tan/5. In the CMSSM, the lowest allowed value for mH+ is found to 
be of the order mfj+ > 400 GeV, while mfj+ ~ 135 GeV is still not excluded in the NUHM 
models. Even lower values could be obtained in the intermediate tan (3 region by relaxing the 
SM bound on rrih. 

It is important not to over-interpret the limits obtained using indirect observables. We have 
shown explicitly that the results from — in particular are subject to large uncertainties 
from the determination of \Vub\- Likewise, the results obtained from K — > yuz/^ are highly 
dependent on the value of fx/ fn from lattice QCD. Improvements in the measurements of B 
physics observables, especially the B^ — > and B DtUt- branching ratios, would certainly 
be welcome and serve to refine the situation. 

We compared the MSSM models, with the constraints applied, to the projected experimental 
sensitivities of ATLAS and CMS in the main charged Higgs discovery channels. This comparison 
illustrates that most of the indirect constraints are relevant in the same parameter space regions 
where the charged Higgs production cross section at the LHC is the largest. We have also 
considered the interpretation of the discovery potential in specific NUHM models. For mH+ < 
rut, we find a possibly large sensitivity to the MSSM benchmark scenario through corrections 
to the bottom Yukawa couphng, while this effect cancels to a large degree for the channel 
r^Uj. when 'mH+ > irit. 

This study can easily be extended to scenarios with other mechanisms than gravity mediated 
supersymmetry breaking. More interesting would be to consider the MSSM beyond minimal 
fiavor violation, including effects of CV- and 7?.-parity violation, or to carry over the constraints 
on charged Higgs bosons to models with enlarged Higgs sectors like the NMSSM [75]. The same 
observables discussed here, in particular those which involve tree-level exchange of the charged 
Higgs, play an important role in constraining any MSSM-like model. We therefore propose 
that complementary discovery channels, governed by couplings which are not constrained at 
this point, should be investigated to determine the prospects for early charged Higgs discovery 
at the LHC. Such a discovery would serve as an indication of a non-minimal model being 
realized in nature. 
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